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In this paper we compute the leading correction to the bipartite entanglement entropy at 
large sub-system size, in integrable quantum field theories with diagonal scattering matrices. We 
O I find a remarkably universal result, depending only on the particle spectrum of the theory and 

not on the details of the scattering matrix. We employ the "replica trick" whereby the entropy is 

■ obtained as the derivative with respect to n of the trace of the n^^ power of the reduced density 
matrix of the sub-system, evaluated at n = 1. The main novelty of our work is the introduction 
of a particular type of twist fields in quantum field theory that are naturally related to branch 
points in an n-sheeted Riemann surface. Their two-point function directly gives the scaling limit 

, of the trace of the n^^ power of the reduced density matrix. Taking advantage of integrability, 

D I we use the expansion of this two-point function in terms of form factors of the twist fields, 

i ' in order to evaluate it at large distances in the two-particle approximation. Although this is 

a well-known technique, the new geometry of the problem implies a modification of the form 
factor equations satisfied by standard local fields of integrable quantum field theory. We derive 

■ the new form factor equations and provide solutions, which we specialize both to the Ising and 
OO . sinh-Gordon models. 
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1 Introduction 



Quantum field theory (QFT) has proven to be one of the most successful theories of the phys- 
ical world. Its main objects are correlation functions of local fields: they describe quantum 
correlations between separated local observables and provide all physical information that can 
be extracted from a model of QFT. In general, correlation functions can be computed only 
perturbatively, giving expressions that apply only in the large-energy region. In contrast, some 
two-dimensional QFTs (one space and one time dimension) have the property of integrability, 
which means that correlation functions are accessible non-perturbatively. In particular, their 
low-energy or large-distance behavior is accessible in exact form in many cases, from exact ex- 
pressions for so-called form factors of local fields [U [2]. This fact has triggered an enormous 
amount of work in computing form factors in a multitude of models of integrable QFT (IQFT) , 
some of which has found applications to low-dimensional condensed matter systems [3] . 

The main characteristic of a quantum system, as opposed to a classical one, is the existence of 
entanglement: performing a local measurement may instantaneously affect local measurements 
far away. This property is essential to the field of quantum computation and teleportation. At 
the theoretical level, there has been considerable interest in formulating measures of quantum 
entanglement [l]-[8] and applying them to extended quantum systems with many degrees of 
freedom, such as quantum spin chains [5]-|17j. One of these measures is entanglement entropy 
[1]. Consider a quantum system, with Hilbert space TL = TLa (^'Hb, in a pure state |^). The 
bipartite entanglement entropy Sa is the von Neumann entropy associated to the reduced density 
matrix of the subsystem A, defined as 

p^ = T¥„,(|V')(Vl), (1.1) 

SA = -TTHApA^Og{pA)) ■ (1.2) 

A recent application of QFT has been to the calculation of entanglement entropy for the 
case of one-dimensional systems |18t 119] . extending primarily on the work [20J, where TIa is 
spanned by the degrees of freedom in some interval A (or set of intervals) of the real line, and 
B is its complement. The authors evaluated the bipartite entanglement entropy in quantum 
systems at criticality, using techniques of conformal field theory (CFT), as well as the leading 
large-distance limit in a massive QFT. For example, when is the ground state and A is an 
interval of length r in an infinite system, they found 

{-log(r/e) e ^ r ^ 
^ 1 (^-^^ 

— -log(emj r 3> m 

where is a correlation length of the QFT, e is some short-distance cutoff, and c is the central 
charge of the CFT. 

In this paper, we will develop a framework for the computation of entanglement entropy in 
massive IQFT using factorized scattering techniques. Our main result is the form of the first sub- 
leading corrections to the entanglement entropy at large distances. Remarkably, we show that 
the leading r-dependent correction is independent of the precise details of the 5-matrix, being 
entirely determined by the spectrum of masses of the IQFT whenever the scattering between 
particles does not involve backscattering. That is, 

1 ^ 

Sa = -\ log(emi) + U - Ko{2rm^) + O (e-^-'-i) (1.4) 
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where a labels the particle types, ma are the associated masses (with m\ the smallest one), V 
is a model-dependent constant, and Kq{z) is the modified Bessel function. For free theories, 
this result was previously obtained by a different approach |2H I22j . The constant U depends 
on the definition of e, but it can be fixed by requiring that no constant correction terms occur, 
for instance, in the upper expression in (jl.3p . With such a definition, we evaluated \J in the 
quantum Ising model, with the result: 

_1 dt f t cosh t 1 £2 

t/ising- glog2-J^ -l^—-^-—^^- — 

which, as we will explain in section O is in very good agreement with existing numerical results 
[llj . This constant could in principle also be recovered by numerically integrating the Painleve 
V equation according to the results obtained in [21j . 

The main novelty of our approach is the introduction of a certain type of twist field T in 
IQFT, whose two-point function is directly related to the entanglement entropy. The initial 
idea [H] is to evaluate the entanglement entropy by a "replica trick" from the partition function 
on a multi-sheeted Riemann surface. The field we introduce naturally arises as a local field in 
an n-copy version of a given model of IQFT, and implements branch points so that its corre- 
lation functions are partition functions on multi-sheeted Riemann surfaces. The large-distance 
corrections to the bipartite entanglement entropy is derived from a large-distance expansion of 
the two-point function of T, obtained by evaluating its form factors. Due to the new geometry 
of the problem, the form factor equations differ from their usual form derived thirty years ago 
[U [2] . The operation of evaluating the entanglement entropy from the result of the form factor 
expansion involves a subtle analytic continuation in n. We provide a general derivation of the 
form factor equations and of the entanglement entropy in the context of diagonal factorised 
scattering theory, and we have studied the Ising and sinh-Gordon cases in detail. 

The paper is organised as follows: in section [21 we discuss partition functions on multi- 
sheeted spaces and their relation to entanglement entropy in general terms. In section [3l we 
develop the form factor program for twist fields, and provide a general formula for the two- 
particle form factors of theories with a single particle spectrum and no bound states, which 
we specialise to both the Ising and the sinh-Gordon cases. In section H] we check our previous 
formulae for consistency by computing the ultraviolet conformal dimension of the twist field 
in the two-particle approximation both for the Ising and sinh-Gordon model. In section we 
compute the two-point function of the twist field in the two-particle approximation and derive 
from that our general formula for the bipartite entanglement entropy for theories with a single 
particle spectrum and no bound states. In section [6] we extend the previous results to diagonal 
theories with many particles and bound states. In section [7| we present our conclusions and 
point out some open problems. Finally we provide four appendices: in appendix |A] we give an 
alternative derivation of the form factors of the twist field in the n-copy sinh-Gordon theory 
using the method of angular quantization; in appendix [B] we compute the vacuum expectation 
value of the twist field in the n-copy Ising model and derive the value (jl.Sp from it; in appendix [Cl 
we present the details of the analytic continuation necessary for the computation of the entropy 
and in appendix [D] we provide the n — > oo limit of the two-point function of the twist field in 
the sinh-Gordon model in the two-particle and saddle-point approximation. 



j = -0.131984... (1.5) 
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Fi gure 1: [Color online] A representation of tlie Riemann surface J^3^ai,i 



2 Partition functions on mult i- sheeted spaces and entanglement 
entropy 

2.1 Partition functions in QFT on multi-sheeted spaces 

The partition function of a model of two-dimensional QFT with local lagrangian density 
£[(/?] (x, y) on a (euclidean-signature) Riemann surface TZ is formally obtained by the path integral 



where [df/^jT^ is an infinite measure on the set of configurations of some field cp living on the 
Riemann surface TZ and on which the lagrangian density depends in a local way. Consider 
Riemann surfaces with curvature zero everywhere except at a finite number of points. Since the 
lagrangian density does not depend explicitly on the Riemann surface as a consequence of its 
locality, it is expected that this partition function can be expressed as an object calculated from 
a model on M^, where the structure of the Riemann surface is implemented through appropriate 
boundary conditions around the points with non-zero curvature. Consider for instance the 
simple Riemann surface M.n,ai,a2 composed of n sheets sequencially joined to each other on the 
segment x £ [ai, 02], y = (see Fig. [T] representing the case n = 3). We would expect that the 
associated partition function involves certain "fields"B at (x, y) = (ai,0) and (x, y) = (02,0). 

The expression (|2.1|) for the partition function essentially defines these fields (that is, it gives 
their correlation functions, up to a normalisation independent of their positions). But in the 
model on M?, this definition makes them non-local. Locality of a field (used here in its most 
fundamental sense) means that as an observable in the quantum theory, it is quantum mechani- 
cally independent of the energy density at space-like distances. In the associated euclidean field 
theory, this means that correlation functions involving this field and the energy density are, as 
functions of the position of the energy density, defined on M? (and smooth except at the positions 
of the fields). The energy density is simply obtained from the lagrangian density, hence it is 
clear that fields defined by (j2.ip in the model on M? are not local. Locality is at the basis of 
most of the results in integrable QFT, so it is important to recover it. 

*Here, the term "field" is taken in its most general QFT sense: it is an object of which correlation functions - 
multi-linear maps - can be evaluated, and which depends on a position in space - parameters x, y that transform 
like coordinates under translation symmetries. 




(2.1) 
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The idea is simply to consider a larger model: a model formed by n independent copies of 
the original model, where n is the number of Riemann sheets necessary to describe the Riemann 
surface by coordinates on M^. Let us take again the simple example of M.n,ai,a2- We re-write 
(f2J]l as 



Mn,ai,a2] = I Wl ' ' ' d^Pnh^ exp 

J C(a-i -.ao) 



dx(iy(£[v9i](x,y) + . . . + £[v9j(x, y)) 



(2.2) 



'C{ai;a2) 

where C(ai, 02) are conditions on the fields ipi, . . . ,ip^ restricting the path integral: 

C{ai;a2) : (/^^(x, 0+) = (/?j+i(x, 0") , x G [ai, 02], « = 1, . . . , n (2.3) 

where we identify n + i = i. The lagrangian density of the multi-copy model is 

£ W [ip^,..., (x, y) = £ [(/^i] (x, y) + . . . + C[ipJ (x, y) 

so that the energy density in that model is the sum of the energy densities of the n individual 
copies. Hence the expression (12. 2 p does indeed define local fields at (ai,0) and (02, 0) in the 
multi-copy model, since this sum is the same on both sides of the segment x E [01,02], y = 
according to the conditions C (01,02). 

The local fields defined in ()2.2p are examples of "twist fields". Twist fields exist in a QFT 
model whenever there is a global internal symmetry a (a symmetry that acts the same way 
everywhere in space, and that does not change the positions of fields): J^2 dxdy C[aip]{x,y) = 
/jg2 dxdj C[ip](x,y). Their correlation functions can be formally defined through the path inte- 
gral: 



(7;(o, b) ■ ■ ■)cM^ oc / [df]u2 exp 



dxdy C[ip](x,j) 



(2.4) 



where • • • represent insertions of other local fields at different positions and the path integral 
conditions are 

Ca-{a,b) : ip{'x.,b~^) = aip{x,b~) , xG [0,00) . (2-5) 

The proportionality constant is an infinite constant that is independent of the position (o, b) and 
of those of the other local fields inserted. The fact that cr is a symmetry ensures that 7^ is local. 
Also, it insures that the result is in fact independent of the shape of the cut in the conditions 
C(j, up to symmetry transformations of the other local fields inserted. A consequence of this 
definition is that correlation functions (7^(o, 6)0(x, y) • • ■)c,r'2 with some local fields 0(x, y) are 
defined, as functions of x,y (smooth except at positions of other local fields), on a multi-sheeted 
covering of with a branch point at (o, 6), whenever aO 7^ O. They have the property that a 
clockwise turn around (o, b) is equivalent to the replacement O aO in the correlation function. 
If aO / O, then O is said to be "semi- local" with respect to 7^. This property, along with the 
condition that 7^ has the lowest scaling dimension and be invariant under all symmetries of the 
model that commute with a (that is, that it be a primary field in the language of conformal field 
theory), is expected to uniquely fix the field 7^, and constitute a more fundamental definition 
than the path integral above, as it does not require the existence of a lagrangian density. We 
will take this point of view in the following, but we will continue to denote a model of QFT by 
C and its n-copy tensor product by 

In the model with lagrangian £("), there is a symmetry under exchange of the copies. The 
twist fields defined by ()2.2p , which we call branch-point twist fields, are twist fields associated to 
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Figure 2: [Color online] The effect of T on other local fields. 



the two opposite cyclic permutation symmetries i i— > i + 1 and i + 1 ^ i {i = 1, . . . ,n, ?i + 1 = 1). 
We will denote them simply by T and T, respectively: 

T = T(j , a : i + 1 mod n 

T = 7^-1 , a^^ : i + i modn 

(see Fig. [5] for the case T) . More precisely, we have 

Z[C, Mn,a,,a^] oc (T(ai , 0)t(a2 , 0)) . (2.6) 

This can be seen by observing that for x G [ai,a2]) consecutive copies are connected through 
y = due to the presence of T(ai,0), whereas for x > 02, copies are connected to themselves 
through y = because the conditions arising from the definition of T(ai, 0) and T(a2, 0) cancel 
each other. 

More generally, the identification holds for correlation functions in the model C on Mn,ax,a2i 
this time with an equality sign: 

. , (T(fli,O)r(a2,O)0.(x,y)---)^(„)^K. 

(0(x, y; sheet i) ■ ■ ■)cM^^^,,a, = 7=7 -r;Er — (2.7) 

(T(ai,0)r(a2,0))^(„)K2 

where Oi is the field in the model iZ^"^ coming from the i^^ copy of C 

The generalisation to Riemann surfaces with more branch points is straightforward, but will 
not be needed here. 

The conformal dimension of branch-point twist fields was calculatecil] in [18]. Consider the 
model £ to be a conformal field theory (CFT). Then also is a CFT. There are n fields Tj{z) 
in vC*-") that correspond to the stress-energy tensors of the n copies of £, and in particular the 
sum T^^\z) = X]j=i '^ji^) stress-energy tensor of C^^\ The central charge of is nc, 

if c is that of C. 

Consider the stress-energy tensor T[w) in C. We can evaluate the one-point function 
{T{w))c,M„ a a by making a conformal transformation from z in R'^ to w in A4n,ai,a2 (here 
z and w are complex coordinates, with for instance z = x + iy) given by 

w — ai 
w - 02 



^In fact, in the paper jlSj branch-point twist fields in the multi-copy model were not introduced explicitly, as 
they are not essential for the evaluation of partition functions in CFT. Only the non-local fields discussed above 
were alluded to, but the method to evaluate the scaling dimension is the same. 
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We have 

where the Schwarzian derivative is 



//\2 



z"'z' -{3/2 ) jz") 
W 



{z,w} 

Using (^(z))^ ]g2 = 0, we obtain 

,\\ - - 1) (ai - a2? 

{-L [W) C,Mn,a-, — n ^ 2 (' \2 / \2 ' 

' '-' ^ 24n^ [w — aij'^lw — a2) 

Since, by (|2.7p . this is equal to {T{ai,0)T{a2,0)Tj{w))^{n)^^2/{T{ai,0)T{a2,0))^(„)^^2 for ah j, 
we can evaluate the correlation function involving the stress-energy tensor of C^^^ by multiplying 
by n: 

(T(ai,0)t(a2,0)r(»)H)^(„)^K. _ c{n^ _ i) (g, _ g^f 

(T(ai,0)T(a2,0))£(„) ]5;2 24n (u; - ai)2(u; - 02)^ ' 

From the usual CFT formula for insertion of a stress-energy tensor 

(T(ai,0)t(a2,0)T(")H)^(„)^R2 = 

Id hi Id h2 \ 

+ 77, — 7^ + — + 77, — 'Tyi W ("i' 0)'^(«2, 0))£(„) 02 



w — a\ da\ {w — ai)^ w — a2 da2 {w — 02)^ 

we identify the scaling dimension of the primary fields T and T (they have the same scaling 
dimension) using (T(ai, 0)T(a2, 0))£(„) ^2 = |ai —02!"^'^": 

dn = ^in-'^ . (2.8) 



12 V 

It may happen that many fields with the main property of branch-point twist fields exist, with 
different dimensions. However, the dimension p.Sp should be the lowest possible dimension. 
Hence, a field with the main properties of branch-point twist field, with this dimension, and 
invariant under all symmetries of the theory should be unique. 

2.2 Entanglement entropy 

Partition functions on Riemann surfaces with branch points can be used in order to evaluate 
the entanglement entropy; this works when A consists of one and also of more than one interval, 
as was explained in [.18j. Consider a (finite or infinite) one-dimensional quantum system with 
Hilbert space TL = TLa ®'Hb where TLa is the space of local degrees of freedom in some interval 
(or set of intervals) A. Consider also the ground state of this quantum system. The bipartite 
entanglement entropy Sa is defined as follows. We first define the induced density matrix as 

PA = Tv-H,m{ij\) (2.9) 

and then we calculate the von Neumann entropy associated to this density matrix: 

SA = -TTHAPA^Og{pA)) . (2.10) 
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This has the interpretation of counting the number of fully entangled "links" between the regions 
A and A as encoded into the ground state {tp). In agreement with this interpretation, it has 
the symmetry property Sa_ = S^, and in a quantum model with local interaction, it is expected 
to saturate to a finite value when both regions A and A are much larger than the correlation 
length. 

The main idea in order to evaluate the entanglement entropy in the scaling limit of a quantum 
model (here, we will consider models on infinite space only) is to use the "replica trick" [18j . 
That is, we evaluate 

T^HaPa (2.11) 

then take the limit n — > 1 of the derivative with respect to n, using the identity log = 
lim„^i l^p^: 

Sa = - lim ^Trn^pl . (2.12) 
n^i an 

This formula requires that we evaluate the trace with real positive n, but the trick is to evaluate 
TrT-^^p^ with positive integer n, then to take the appropriate "analytic continuation" (it is 
unique if we assume a certain asymptotic behaviour as n ^ oo - we will discuss this in section 
[5]). Considering positive integer n is useful, because in the scaling limit, denoting the QFT 
model associated to the region near the critical point by £, we have, taking A to consist of only 
one interval, 

TTn^pl^Z[C;Mn,a„a,] (2.13) 

where oi and 02 are the dimensionful end-points of the region A (the scaling limit is taken with 
the length of the region A in proportion to the correlation length, which is then sent to infinity). 
As we saw above, this can be computed as a two-point correlation function of local fields in 
£(") using (12. 6p . More precisely, with m a mass scale of the QFT and e the some dimensionful 
distance of the order of the site spacing, me being in inverse proportion to the dimensionless 
correlation length, we have 

Tin.pl ~ ^ne2^"(T(ai,0)t(a2,0))^(„)^K2 (2-14) 

with an n-dependent non-universal normalisation constant Zn (with Zi = 1), and where dn is 
the scaling dimension (j2.8p . For later convenience, e is chosen in such a way that dZn/dn = 1 
at n = 1. Note that the expression above is dimensionless, since the operators T and T both 
have dimension dn (in particular, their individual vacuum expectation value in the QFT is 
proportional to m'^" ). 

Similarly, for regions composed of many disconnected components, Tr 7.^^/3^ is identified with 
partition functions on Riemann surfaces with many branch points, as explained in [18], but we 
will not consider this case here. 

3 The form factor program for branch-point twist fields 

We now turn to the description of QFT on Minkowski space-time in terms of its Hilbert space 
of asymptotic relativistic particles. In the context of 1 -|- 1-dimensional IQFT, form factors are 
defined as tensor valued functions representing matrix elements of some local operator 0{x) 
located at the origin x = between a multi-particle in-state and the vacuum: 

F^^^-^^e,, ...,6k) := mme,, . . .,e,y;^_^^ . (3.1) 

Here |0) represents the vacuum state and \9i, . . . , OkY^^^,,,^^^ the physical "in" asymptotic states 
of massive QFT. They carry indices /ij, which are quantum numbers characterizing the various 



7 



particle species, and depend on the real parameters 9i, which are called rapidities. The form fac- 
tors are defined for all rapidities by analytically continuing from some ordering of the rapidities; 
a fixed ordering provides a complete basis of states. 

The main characteristics of massive integrable models of QFT are that the number of particles 
and their momenta set are conserved under scattering, and that the scattering matrix factorises 
into products of two-particle scattering matrices, which are the solutions of a set of consistency 
equations and analytic properties. These consistency equations and analytic properties are often 
strong enough to completely fix the scattering matrix in integrable models. Similarly, the form 
factors are fixed by a set of equations and analytic properties depending on the two-particle 
scattering matrix (or S-matrix). 

In this section we want to show how the standard form factor equations for 1 + 1-dimensional 
IQFTs must be modified for the branch-point twist fields. Let us consider an integrable model 
consisting of n copies of a known integrable theory possessing a single particle spectrum and no 
bound states (such as the Ising and sinh-Gordon models). We have therefore n particles, which 
we will denote by indices 1, . . . , n. The S'-matrix between particles i and j with rapidities 9i and 
9j will be denoted by Sij{9i — 9j) (that it depends on the rapidity difference is a consequence of 
relativistic invariance). Particles of different copies do not interact with each other, so that the 
S'-matrix of the model will be of the form 

Sii{9) = S{9) V i = l,...,n, (3.2) 
Sij{9) = 1, V i,j = l,...,n and i ^ j, (3.3) 

where S{9) is the 5-matrix of the single-copy integrable QFT. As explained above, as a conse- 
quence of the symmetry of the model, a twist field T must exist such that if ^'i, . . . , are the 
fields associated to the fundamental particles of each copy of the original model, then the equal 
time {x^ = y^) exchange relations between T and ^'i, . . . , ^„ can be written in the following 
fornil 

■^i{y)T{x) = T(x)^'i+i(y) > y\ 

^^{y)T{x) = T{x)^i{y) x' < y\ (3.4) 

for i = 1, . . . ,n and where we identify the indices n + i = i. The relation with the previous 
section is obtained by recalling that in going from the Hilbert space description to the path 
integral description, the order of operators is translated into time-ordering (or y-ordering in 
euclidean space), and that left-most operators are later in time. It is well known that such 
exchange relations play an important role in the derivation of the consistency equations for the 
form factors of the operator T. Generalising the standard arguments to the exchange relation 
()3.4p . the form factor axioms are 



/ |---A'iA'i+i--- fan \ 

^r|MiM2...M.(^^^27ri,...,^fe) 
-iResg^^e^F^f^--^^ (^0 + ivr, 9o,9,..., 9k) 

-iRes-s.F^lf'''---^'^{9o + ivr, 9o,9,..., 9u) 



— Q (Q \Tp'■\^^■^^i+l^^i■■^( O O \ 

— ^l^il^i+iK^ii+lJ^k (.• • • ,t/i+l,t/i, • • -j, 

— 1^2, • • • , Pfc, t7ij 

= FP^--^''{9u...,9k), 

k 

= -llS>^^^.((^o^)F^^'■■■^'{0l,■■■,Ok). (3.5) 

i=l 



■■"Here we employ the standard notation in Minkowski space-time: a;" with u = 0,1, with x being the time 
coordinate and being the position coordinate. 
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Figure 3: [Color online] A pictorial representation of the effect of adding 27rz to rapidity 6i in 
form factors of the twist field T. 



Here 9ij = 6i — 9j and the first axiom is in fact the same as for local fields. In the second 
equation, the crossing or locality relation, we introduced the symbols /xj = //j + 1. As compared 
to the usual form factor equations, it is altered by the nature of the exchange relation and it 
now relates form factors associated to different particle sets. Finally, the last two equations 
generalise the standard kinematic residue equation to branch-point twist fields. Once more, 
the exchange relations (13. 4p are responsible for the splitting into two equations. Here, for later 
convenience, we wrote the equations in their general form valid also for many-particle models, 
where Jl represents the anti-particle associated to /x. In the present case, the integrable model 
we started with has just one particle (so that /i labels the copies) and therefore each particle is 
its own anti-particle. Since we are dealing with theories without bound states, these are in fact 
all the equations which one needs to solve. Pictorial explanations of the second and of the last 
two equations are given, respectively, in Figs. [3]and[ll 

The other field T introduced in section [2] is a twist field with similar properties as T but 
whose exchange relations with the fundamental fields of the theory are given by 

^^{y)f{x) = f{x)^^{y) x^<y\ (3.6) 

instead of (j3.4p . This implies that on the Hilbert space, we have 

f = TK (3.7) 

In order to fully define the fields T and T, we need to fix their normalisation, which does 
not follow from the form factor equations. We will adopt the usual CFT normalisation: 

(T(x)t(O)) ~ r-^"'" as r ^ . (3.8) 

Here and below, r denotes the space-like separation (x^)^ — (a;^)^, and the two-point function 
just depends on it thanks to relativistic invariance and spinless-ness of the fields involved. 

3.1 Two-particle form factors 

As usual in this context, we define the minimal form factors f'^I^'^ , n) to be solutions of the 
first two equations in (j3.5p for k = 2 without poles in the physical sheet lm(9) S [0, tt]. That is, 

= Fl\i\-e,n)S,,{e) = F^i'^\2ni -e,n) V (3.9) 
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Figure 4: [Color online] The kinematic poles come from the structure of the wave function 
far from the local fields, at positive and negative infinity. Adding in to rapidity 6i puts the 
particle in the "out" region. With a particle in that region, there are delta-functions representing 
particles in the "in" region going through without interacting with the local fields. Those occur 
from the e*^^ form of the wave function at positive and negative infinity. But if the coefficients at 
both limits are different, S-e^^^ and S+e^^^ with S- ^ S+, then there are also poles in addition 
to these delta- functions. Only these poles are seen in the analytic continuation 9i ^ 9i + m. 
Different coefficients come from the semi-locality of the twist field and the non-free scattering 
matrix, as represented here. 



where the S'-matrix is given by (|3.2p . Repeated use of the above equations leads to the following 
constraints: 

FT^'iO^n) = F2^^\9,n) V ^,J,k (3.10) 
^Iin'(^>^) = F^S.\27Tij-l)i-e,n) V j/1. (3.11) 

"Till 

These equations show that computing just the form factor F^^^ {6, n) is enough to determine all 
minimal form factors of the theory. A consequence of these equations is that this minimal form 
factor must have no poles in the extended strip lm(0) G [0, 27rn]. From the equations above it 
is easy to deduce 

^nTin'(^>^) = Fl!^{-(^^ri)S{e) = Fl}^{-e + 2^m,n). (3.12) 

In order to develop a systematic procedure to solve these equations it is useful to recall that, for 
a standard local operator the minimal form factor equations take the form 

hi{e) = hi{-0)S{ne) = fui-9 + 27ri), (3.13) 

provided that the S'-matrix of the theory is given by S{nO). Thus given a solution to the previous 

'7~l 11 

equation, the function F^^^ {0,^) = fii{9/n) is automatically a solution of ()3.12p . 

In the context of integrable models, a systematic way of solving such type of equations has 
been developed whereby, given an integral representation for S{6), an integral representation 
of fii{9) can be readily obtained [1]. For diagonal theories, the integral representation of the 
S'-matrix takes the form 

S{9) = exp / —g{t) smh - 



ITT 



(3.14) 
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7~l 11 

Figure 5: [Color online] The structure of the function (^) in the extended sheet Im(^) G 

[0, 27rn], in the case n = 3. Crosses indicate the positions of the kinematic singularities. Shaded 
regions represent the physical sheets of the form factors F^^^'^iQ^^ F^^'^'^iQ) and F^^^^{6). 



where g{9) is a function which depends of the theory under consideration. A trivial consequence 
of the previous equation is 



S{n6) = exp 
and from here, it is easy to show that 
hi{9) =AAexp 



f dt , , ^ . rt9 
Jo t V^^ 



dt 



-g{t) sin^ 



itn 



1 + 



/g tsinh(nt) ' y 2 y vr 
where is a normalization constant. Therefore, the desired solution is 



,n) = hi{e/n) = AAexp 



dt 



-g{t) sin^ 



n + 



vr 



(3.15) 



(3.16) 



(3.17) 



/o tsinh(nt)' 

So far in this section we have computed the minimal form factors. However what we ultimately 
need are the full two-particle form factors. These are solutions of ()3.9p which include poles in the 
extended physical strip mentioned before. Their pole structure is determined by the kinematical 
residue equations together with (13. 9p . According to these equations the form factor F2 {6,n) 
has two poles in the extended physical sheet at 9 = itt and 9 = i7r(2n — 1). It is not difficult to 
show that a solution of (13. 9p which is consistent with the above pole structure is given by 



(0) 



(r)sin^ 



FZ 



\9,n) 



2n sinh 



iTr{2{j-k)-l)+e 
2n 



sinh 



»7r(2(fc-j)-l)-e 
2n 



F. 



r\jk 



(ivr, n) 



(3.18) 



where the normalization has been chosen so that the kinematical residue equation gives 

(T). (3.19) 



In addition, the constant factor sinh(7r/n) guarantees that all form factors vanish for n = 1 as 
expected, since in that case the field T can be identified with the identity. The structure of 
two-particle form factors is depicted in Fig. [5l 

For the field T, the exchange relations imply that form factors of the field T are equal to 
those of the field T up to the transformation i ^ n — i for each particle i. At the level of the 
two particle form factors, this means that 



FP^{9,n) = F, 



f\{n-i){n-j) I 



,n 



(3.20) 
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This property can be combined with (j3.10p - (j3.1ip to show that 



F^'''{e,n) = ^2 (^,n). (3.21) 
Fp^^{e,n) = Fp^\e + 2TTi{j -l),n). (3.22) 



3.2 The Ising model 



The Ising model is, together with the free Boson theory, the simplest integrable model we can 
possibly consider. It describes a real free fermion and therefore the scattering matrix is simply 

S{e) = -1. (3.23) 

Form factors of local fields of the Ising model were first computed in [U |23] and later on in [24j 
for so-called descendant fields. A solution of (j3.9|) for j = A; = 1 is given by 

^2'W = -^smh(|-). (3.24) 

This is in fact the standard minimal form factor already employed in [H [23], with 9 — > 9/n. 



3.3 The sinh-Gordon model 



The sinh-Gordon model is a quantum integrable model possessing a single particle spectrum 
and no bound states. The corresponding S'-matrix [25l [261 12Z] is given by 



s{e) 



tanh ^( 



^^) 



tanh ^ ( 



(3.25) 



The parameter B € [0,2] is the effective coupling constant which is related to the coupling 
constant (3 in the sinh-Gordon Lagrangian [281 [29j as 



2/32 



(3.26) 



87r + ' 

under CFT normalization [30]. The S-matrix is obviously invariant under the transformation 
B 2— i?, a symmetry which is also referred to as week-strong coupling duality, as it corresponds 
to B{P) B{8ttP^^) in (I3.26p . The point i? = 1 is known as the self-dual point. Form factors 
of the sinh-Gordon model were first computed in |3T]. The program was thereafter extended 
to other operators in [32j and more recently in [33j. The S'-matrix above admits an integral 
representation which is given by (|3.14p . with 



8sinh(^f)sinh(| (l 



9it) 



Therefore, the minimal form factor is given by 



^))sinh(i) 



sinht 



^Ii'n'W=exp 



dt sinh ^ sinh . / it 

- cosh t [ n + 



t sinh(nt) cosh 



vr 



(3.27) 



(3.28) 
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where we have chosen the normahzation J\f 
dt smh(at) sinh(/3t)e~'^* 



7"l 11 

-^min [i'^iT')- Employing the identity 



t 



sinh(tit) 



- log 

2 ^ 



a+/3+7+M 
2u 



-g— /3+7+M 
2u 



2u 



2m 



(3.29) 



where T{x) is Euler's gamma function, we obtain the alternative representation 

^ 2n-2w+B+2k \ -p f 2n+2w±B+2k'' 



oo 

E 

fc=0 



(-l)'^log 



4n 



4n 



9 



■p / n—w+k 
^ \ 2n 



271 



■p / 2n-2^>+2-_B+2fc ^ p / 2n+2w+2-_B+2fc ^ 



n— ui+fc+l 
2n 



n+ui+fc+l ^ 
2n ) 



(3.30) 



with w = n + iO /it, OT equivalently 



n 

fc=0 



' 2n-2M;+_B+4fc 

4?! 



)r( 



2n+2w+B+Ak 
in 



)r( 



2n-2M;+2-_B+4A; ^ p ^ 



2n+2M;+2-_B+4fc ^ 
4n ) 



p / n— ui+2fc 
I 2n 



n+^>+2fc ^ 
2ri ; 



p / n-ui+2fc+2 
I 2n 



)r( 



n+^)+2fc+2 ^ 
2n y 



' 2n-2«;+-B+4fc+2 ^ 
4n ; 



p ^ 2n+2w+_B+4fc+2 ^ p ^ 2ra-2w+4-_B+4fc -j p ^ 



2n+2^>+4~_B+4fc ^ 
4n ; 



(3.31) 



As a consistency check, it is quite easy to show that for n = 1 the minimal form factor above is 
the standard minimal form factor associated to local fields in the sinh-Gordon model computed 
in [31]. In appendix A we will show how the same expression can be derived form the angular 
quantization scheme proposed in and later carried out for the exponential fields of various 
models (including the sinh-Gordon model) in j35j. 



4 Identifying the ultraviolet conformal dimension of T 

In this section we verify that the form factors constructed above agree with the properties of the 
operator T at conformal level, that is, in the ultraviolet limit. As is well-known, the form factor 
program provides a way of carrying out this verification by allowing us to compute (at least in 
an approximate way) the correlation functions of various fields of an integrable quantum field 
theory. In the ultraviolet limit, it is possible to relate a particular correlation function to the 
holomorphic conformal dimension of a primary field by means of the so-called A-sum rule: 

^ ^ ^e{r)f{0))dr (4.1) 



A 



r 



2(T) ,0 

(where the integration is on a space-like ray), originally proposed by G. Delfino, P. Simonetti and 
J.L. Cardy in [36], where G is the local operator corresponding to the trace of the stress-energy 
tensor. The first equality, expected from CFT, holds from the A-sum rule thanks to the fact 
that B commute with T and that 0^ = B. The holomorphic conformal dimension is related to 
the scaling dimension by dn = 2A^ , where dn is expected to be (|2.8p . 

By introducing a sum over all quantum states and carrying out the r-integration, the ex- 
pression above can be rewritten as 



1 



2(T) 



oo 

E E 

k=i Mi-.-Mfc. 



d9i... dOi, 



0|Mi---Mfe/ 



c) F> 



fc!(2^)ME: 



m^. cosh 6, 



(4.2) 
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where the sum in /ij with i = 1, . . . ,k is a sum over particle types in the theory under consid- 
eration. The sum starts at A; = 1 since we are considering "connected" correlation functions, 
that is, the A; = contribution has been subtracted. The sum above, can only be carried out in 
particularly simple cases. For most models, one must be content with evaluating just the first 
few contributions to the sum. Fortunately, the many studies carried out in the last years provide 
strong evidence that the sum above is convergent and that in fact, the first few terms provide the 
main contribution to the final result. Indeed, the convergence is often so good that considering 
only the contribution with k = 2 already provides very precise results (see e.g. |31|). Expecting 
a similar behaviour also in our case, we will approximate the sum above by the two-particle 
contribution, that is 



2{T) J J 2(2^)2m2 (cosh 01 + cosh ^2)^ ■ 



The factor of n is a consequence of summing over all particle types and using (j3.10p . In addition, 
the only non-vanishing contribution comes from form factors involving only one particle type, 
since we are considering n non-interacting copies of the model. This implies that 

F^\'\e) = Q V i/j. (4.4) 

Changing variables to 6 = 9i — O2 and 0' = 61 + 62 we obtain, 

dede'F®\^^{e)F'^^^^{e,ny 



n 



'^iV J J 2(2^)2m2 (2cosh(0/2)cosh(072)) 



2 



-00 — CXD 

00 



n 



F?'\e)Fr'\e,nr 

327r2m2(T) J cosh2(0/2) ' ^ ' 

—00 

Let us now evaluate this integral both for the Ising and sinh-Gordon models. 
4.1 The Ising model 

For the Ising model the only non- vanishing form factor of the trace of the stress-energy tensor 
is the 2-particle form factor. Hence, the two-particle approximation (|4.3p becomes exact. The 
two-particle form factors are given by 

fP'\0) = -i(T) cos(^) / ^\ eiii(^) ^ _ 2 ^.^^ ( 0\ ) 

2 ^' ^sinh(^)sinh(^) \2nJ' ^ ^' \2j' ^ ' 

and therefore 

cos (^) sinh il) sinh (f ) 



167r y sinh (^) sinh (^) cosh 



2 /e- 



■de. (4.7) 



It is easy to check numerically that the above integral exactly reproduces the expected value 
for c = 1/2, 

2A- = l(n-l)=d„, (4.8) 
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for any value of n. The integral can also be computed analytically, at least for n even. In this 
case, shifting t by 27rnz the integral above changes by a sign and therefore it is possible to show 



2A^ = 27ri ^ : 



J' 



(4.9) 



where rj are the residues of the poles of the integrand at t = ivr(2j — 1), with j = 1, . . . ,n. Of 
those, the poles at j = l,n are triple whereas all the others are double poles. A tedious but 
straightforward computation yields 



27riri 



27rzr„ 



-l + n2 + 6cot(^)/sin(^) 



48n 

(cot((^)/sin(i^) +cot(f )/sin(f ) 
8n 

(_l)n+i (4 _ ^2 _ 3cot2(^) + 6cot^(^)) 
48n ■ 



1 <j <n, 



(4.10) 



Finally, we need to add up all these residues. The sum over the rj residues becomes in fact very 
simple, since it is a telescopic series. We obtain. 



n-l 

27ri rj 
i=2 



cot( 



4nsin(^)^ 



(4.11) 



which gives 
4.2 The sinh-Gordon model 



In this case, the relevant 2-particle form factors are given by 

nTIll, 



Fr 



rill 



(0) 



(T) sinh I 



2n sinh 



2n 



sinh I 



,T|11 
min ^ 



Fr 



eiii 



2Trm^I^^-i^. (4.12) 



The form factors of were computed in [3T]. Since is a local operator, its minimal form 
factor is given by (|3.3ip with n = 1. Thus, 



A 



T 



sm 



F2'{0,ny F^( 



1 

I -nh(^)sinh(^)cosh^(f)F2\.vr,n)*F2\.vr,l) 



dO. 



(4.13) 



The tables below show the result of carrying out this integral numerically for various values of 
n and B. Next to each value of n in brackets we show for reference the expected value of A-^, 
as predicted by the CFT formula ()2.8p (with, again, A-^ = d„/2). 





n = 2 (0.0625) 


n = 3 (0.1111) 


n = 4 (0.1563) 


n = 5 (0.2) 


B = 0.02 


0.0620 


0.1114 


0.1567 


0.2007 


B = 0.2 


0.0636 


0.1135 


0.1599 


0.2048 


B = 0.4 


0.0636 


0.1148 


0.1620 


0.2074 


B = 0.6 


0.0643 


0.1155 


0.1631 


0.2088 


B = 0.8 


0.0644 


0.1158 


0.1636 


0.2096 


B = 1 


0.0644 


0.1159 


0.1637 


0.2098 
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n = 6 (0.2431) 


n = 7 (0.2857) 


n = 8 (0.3281) 


n = 9 (0.3704) 


n = 10 (0.4125) 


B = 0.02 


0.2436 


0.2864 


0.3289 


0.3712 


0.4135 


B = 0.2 


0.2488 


0.2925 


0.3360 


0.3793 


0.4225 


B = 0.4 


0.2522 


0.2966 


0.3407 


0.3846 


0.4284 


B = 0.6 


0.2540 


0.2988 


0.3433 


0.3876 


0.4317 


B = 0.8 


0.2550 


0.2999 


0.3446 


0.3890 


0.4334 


B = l 


0.2552 


0.3002 


0.3449 


0.3895 


0.4339 



The figures obtained are extremely close to their expected value for all choices of B and n. In 
most cases they are slightly above the expected value. This is not surprising since the 4- and 
higher particle contributions are not necessarily positive. 



5 Two-point functions and the entanglement entropy 



As explained before, the entanglement entropy is given by the derivative with respect to n of 
the two-point function (T(r)T(O)) evaluated at n = 1. The behaviour of the entropy at short 
separations r <C is described by the conformal limit of the model and is already well known 
[181 119j . At large separations r ^ (in the infrared limit) it is also known to saturate. Here 
we would like to evaluate the first correction to the entropy in the infrared limit. In this limit, 
the two-particle contribution provides the first sub-leading exponential term in the correlation 
function (T(r)T(O)). Hence, the two-particle approximation should provide both the saturation 
value of the entropy, coming from the disconnected part of the correlation function, and the 
exact first exponential correction at large rm, coming from the two-particle contributions. 
The two-point function in the two-particle approximation is given by 



OO — OO 

OO 



^— rm(cosh 6i+cosh 62) 



(r)Ml + ^ J def{e,n)Ko{2rmcosh{e/2)) 



(5.1) 



where we changed variables as in section HI Kq{z) 
out one of the integrals, and we defined 



is the Bessel function resulting from carrying 



(5.2) 



n-l 



n 



(with f{9, 1) = 0). Notice that the function above is only defined for integer values of n. 

In order to obtain the entropy we should now analytically continue the two-point function 
()5.ip . as function of rm and n, from n = 1, 2, 3, . . . to n G [1, oo), compute the derivative with 
respect to n and evaluate the result at n = 1. The analytic continuation is of course not unique. 
We will choose the one which is such that for Re(n) > 0, the two-point function divided by (T)^ 
is 0{e'^"') as n — > oo for some q < n. It is then unique by Carlson's theorem [37]. This choice is 
motivated by the fact that the trace (j2.1ip has this behavior with g < for any finite system. 
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since the eigenvalues of are real an positive, and are normalised to sum to 1. Of course, 
this is merely a motivation. For infinite systems, eigenvalues should have dense components and 
we could have algebraic behaviors; the scaling limit and the limit n — > oo may not commute, 
so that the coefficient of the small-distance power-law e^'^" may be divergent; and the limit 
rm — > oo and n — > oo of this coefficient may also not commute, so that the behavior in n of 
the large-rm expansion coefficients may also be divergent. We expect that this gives at most 
algebraic divergencies in n, but a better understanding would be desirable. 

There are three main observations necessary to understand the analytic continuation of (|5.1|) 
and the evaluation of the derivative at n = 1: 

• Structure of the analytic continuation: There is no natural (as described above) analytic 
continuation of /(O, n) from n = 1, 2, 3, . . . to n G [1, oo]. Instead, there is such an analytic 
continuation from n = 2, 3, 4, . . . to n G [l,oo], which we will denote by f{n), but it has 
the property that 

/(l)y^/(0,l)=0. (5.3) 

On the other hand, the function f{9,n) does not have any such feature for 9^0: its 
analytic continuation f{9,n) to n £ [l,oo] agrees with f{9,n) for all n = 1, 2, 3, . . .. In 
particular, f{9, 1) = 0. Hence, as n — > 1, the function f{9, n) does not converge uniformly 
on G (— oo, oo). 

• Kinematic singularities: The non-uniform convergence of f{9, n) can be seen to be a 
consequence of the collision of the kinematic singularities of form factors {9, n) (at 
9 = in and 9 = i7r(2n— 1)) that occur when n — > 1 (see Fig. [6]). The idea is made clear from 
considering Poisson's resummation formula. Consider a sum of the type X^"=i s{9,j), as in 
the r.h.s. of (15. 2|) . re- written using Poisson's resummation formula, which holds whenever 
s{9,n) = s{9,0) (the summand s{9,j) in (|5.2p indeed satisfies this by unitarity): 

'^s{9,j)= ^(s„fc-Sfc), Sfc = / dje ^s{9,j) . (5.4) 

j=l k=-oo 

For any 7^ 0, no singularity of s{9,j) occur on the integration pathH defining s^, so that 
when n — > 1, the result should vanish. In fact, since s{9,j) itself vanishes like (n — 1)^ 
for 9^0, the result vanishes like (n — 1)^. When — > 0, the poles of s{9,j), which are 
at j = ^ ± 2^ and j = n — ^ ± pinch the integration path, but the quantities Sk 
are still finite because the divergent contributions at Re(j) = ^ and Re(j) = n — | cancel 
out - we are left with principal-value integrals with two double-poles (and the resulting 
conditionally convergent sum over k has a unique finite value defined by the limit 9^0). 
However, since these two double poles collide at n = 1 and fuse into a higher order pole, 
there is no guarantee that the result vanishes as n — > 1. 

• delta-function: From the arguments above, the full contributions of the derivative with 
respect to n of the two-particle contribution to the two-point function, in the limit n ^ 1, 
is obtained from the region ~ 0, and is due to the kinematic singularities. In fact, notice 
that as re — > 1, the derivative with respect to re of f{9,n) "around" 6* = should diverge: 

''The sum over k can be made absolutely convergent by a slight imaginary shift of the j-integration path 
defining Sk- Since s(6,j) is in fact not periodic in j, there are contributions along Re(j) = Q,n additionally to 
the shifted path, with coefficients that vanish linearly as Im(j) 0. This linear vanishing guarantees a vanishing 
like l/k^ of Sk at large fe, making the sum convergent. 



17 



27rin 



>^27rin-7ri 



n^l 



27rin 



,11 



711 
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Figure 6: [Color online] The collision of kinematic singularities when n — > 1. 

indeed, the values of f{9,n) at 9 = and at its neighboring points should be about the 
same for all n > 1, by continuity, but they should reach a finite separation at n = 1; hence 
variations must be very strong near to n = 1. This leads to the expectation that -§^f{9, n) 
at n = 1 is proportional to 6{6). The contribution of the kinematic singularities to the sum 
in the function f{9, n) (j5.2p is obtained from the singular behavior in j of the summand 

s{9,j)=F^\^\-9 + 2mj,n) [f^^^^)* {-9 - 2^ij,n): 

,^ , zF,^'"(-26' + 27rm-m) i ^^'"(-26* + ivr) 

s(9,j) ~ = - = - - + c.c. 

' —9 — 2TTij + 2-Kin — m —9 — 2TTij + in 

(where c.c. means complex conjugate, for real 9). It is a simple matter to perform on this 
expression the sum Yll=i^ giving 



n-l 



:T\U, 

ITT) 



+C.C. 

where ip{z) = dlogT{z)/dz is the derivative of the logarithm of Euler's Gamma function. 
This has no poles at 9 = 0, as the kinematic poles of the form factors involved cancel out. 
The poles that are nearest to Re(^) = as n ^ 1 are at 9 = ibz7r(n — 1), coming from the 
form factors involved. The residues to first order in n — 1 give: 

N ::,-.( inin — 1) zvrfn — 1) \ , , , , 

/ 0, n ~ / I \ - \ \ n ^ 1 5.5 

\2[9 + nr[n — Ijj 2[9 — nr[n — Ijj / 

with 

/(I) = I ■ (5.6) 

This has simple poles at 9 = ibi7r(n — 1) with residues that vanish at n = 1, gives /(I) at 
^ = and vanishes like {n — 1)^ as n ^ 1 for 9^0. The limit n — > 1, as a distribution on 
9, is easily evaluated: 

I^Mn)) _^ = n'f{m9). (5.7) 

In Appendix [Cl we give the full form of f{9,n) and verify that this is correct. Note that 
for the free case our result is in agreement with the n — > 1 limit evaluated in |21j . 
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Inserting this inside (|5.ip gives, using (j2.14p and (j2.12p . the entanglement entropy SA{i"m) 
for A an interval of length r: 

SA{rm) = -- log(em) + U - -Ko{2rm) + O (e~^™) (5.8) 
3 8 

where 

an V /n=i 

Therefore the sub-leading large rm terms in the entropy are given by the Bessel function 
KQ[2rm), up to terms that are exponentially smaller (coming from the neglected 4-particle form 
factors). The term U involving the derivative of the vacuum expectation value of T has a univer- 
sal meaning since the normalisation of T has been fixed (see p.8p ). More precisely, with this nor- 
malisation, the entanglement entropy at short interval length is S'A(rm) = — | log(e/r) + 0(rm). 
It is of course possible to define the quantity [/ in a way that is obviously universal, valid for 
any choice of short-distance cutoff: 

U = hm (Sa{0 - SAiC^) - ^ \ogi) . (5.10) 

We will now proceed to evaluate explicitly the functions f{0,n) and f{n) for the Ising and 
sinh-Gordon models, verifying some of the results above. It is worth noting that the function 
/(O, n), which we study in more detail below, also has a meaning as the coefficient of the leading 
exponential correction to the partition function on the Riemann surface Mnfl,r' 

(™,.,r,(i.=Z<^.o(0)) 

5.1 The Ising model 

For the Ising model, the function ()5.2p at = is given by 



i=2 sm" [^^2^ ) [^^2^ ) 



The result of the sum can be obtained analytically as in appendix O In this case, the integral 
part of the formula ()C.6p is zero, so that f{n) = 1/2 for all n. That is, the connected part 
of the correlation function at large rm behaves linearly with n, for all values of n. Hence, the 
entanglement entropy can be computed to 

SAirm) = log(em) + ^Jismg - ^Ko(2rm) + O (e"^^"^) . (5.13) 

The constant ?7ising can also be evaluated explicitly using the relation between the Ising model 
and the free Dirac fermionic model, as explained in Appendix [Bl We find that the expectation 
value of the branch-point twist field is 



m\ 24 



{T) = (^)-' "exp 



dt f I n e~2* / 1 

n 

n 



4t V sinh t sinh ^ sinh^ t 



n 



(5.14) 
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which gives 



^ lo 2 n (^'^^^^^ 1 
6 °^ ~ io 2t V sinh^ t ~ sinh^ t 



f/ising = - log 2 - / - ( — ^ - — ^ - ^ ) = -0.131984... (5.15) 



Once the value of Uismg has been fixed we can compare our expressions for the entropy in the 
deep-infrared and deep-ultraviolet regimes to existing results in the literature for the quantum 
Ising chain. The comparison goes as follows: from [19ll38j it is possible to obtain the expression 
of the entropy at large separations r S> m^^ in terms of the lattice spacing a 

5A = -^log(am) + ^log2. (5.16) 

By comparing (j5.16p to our formula 

Sa = log(em) + f/ising, (5.17) 

we obtain the precise relationship between the short-distance cutoff e and the lattice spacing a. 
That is 

e = ^e^^i-"s = (0.0566227.. .)a. (5.18) 

8 

As mentioned above, formula (j5.16p follows from the results in [19\ I38j . However this is 
not entirely trivial and some clarifications are due here. First, the formulae given in these 
publications are expressed in terms of a parameter k rather than a. The parameter k is related 
to the value of the transverse magnetic field h of the quantum Ising chain as follows: 

where h = 1 corresponds to the critical point. The entropy in both regions of values of k was 
computed in [38] whereas in |19J only the h < 1 regime was considered. However, it is easy 
to show that both regimes give the same infrared value of the entropy when h approaches 1. 
Second, it is a standard result that itam = 1 — h, where the positive sign corresponds to /i < 1 
and the negative sign corresponds to h > 1. Therefore, the fact that the infrared value of the 
entropy is the same both for h > 1 and /i < 1 is in agreement with what we expect from QFT, 
as going from one regime to the other amounts formally to a change in the sign of the mass 
m which has no effect on the value of the energy density. That is, the two regions should be 
described by the same QFT in the scaling limit. The relations between k and h and between 
am and h given above allow us to relate the parameter k (in terms of which the entropy is 
expressed) to the lattice spacing a. Once this is done we only need to expand the expressions in 
[191 138] around the value k = 1. Since these formulae are given in terms of the complete elliptic 
integral of the first kind K(k), we need to employ the standard expansion 

K{k) l/21og(l - /c) -F3/21og2 for k ^ 1, (5.20) 

to obtain (j5.16p . 

Plugging (I5.18P into the short-distance expression of the entropy we obtain: 

Sa = -- log(e/r) + 0(rm) = -- log(a/r) + 0.478558... + 0(rm) . (5.21) 
6 6 
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Figure 7: [Color online] A plot of the function Sl = \ \og2{L) + 0.690413 for L < 200. The 
graph is in very good agreement with the numerical values plotted in Fig. 9 of [llj. 

We can compare this result to the numerical values obtained in [TT] for the Ising spin chain. 
However, we must notice first that in there the entropy was defined as 

SL = -TTHApA^Og2{PA)) , (5.22) 

whereas in this paper we have used the definition (II. 2p . This means that in order to compare 
(|5.21|) to the results obtained in |rl] we must divide our formula by a log 2 factor. This gives 

Sl = - \ogniL) + 0.690413... (5.23) 
6 

where we now introduced the parameter L = r/a which is the number of sites in the interval 
A. A plot of this function is presented in Fig. [7| which is to be compared to the blue curve in 
Fig. 9 of [11]. Very good agreement is found, which supports the twist-field realization proposed 
in appendix [Bl 

It is worth noticing that in the free case there are alternative ways of computing the entan- 
glement entropy [39l SOl HI] . 

5.2 The sinh-Gordon model 

In this section we will verify that /(I) = 1/2 also for the sinh-Gordon model. Here, the integral 
part of (jC.6p is not vanishing, hence in contrast to the Ising model result (|5.12p . f{n) is in 
general not constant with n. This can be easily seen from Fig. [8] where the functions nf{0,n) 
and n/(oo) are presented, both for the Ising and sinh-Gordon models. 

The non-linearity of the function nf{n) in the sinh-Gordon model can be seen more clearly 
by magnifying the lower left corner of Fig. [8] (a), which we have done in Fig. [8] (b). Figure (b) 
also shows how /(I) appears to tend to the value 1/2 and how the function nf{n) deviates from 
the straight line for small values of n. Both for the Ising and sinh-Gordon model the figures 
show clearly that the function n/(0,n) has a jump at n = 1, since /(0, 1) = 0. We will now 
attempt to provide a more rigorous description of these behaviours. 
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Figure 8: [Color online] (a) shows 4 functions: the points are the function n/(0,n) for integer 
values of n in the interval [2, 50] both for the Ising (black squares) and sinh-Gordon (black circles) 
models, evaluated numerically. The solid blue line gives the corresponding analytic continuation 
nf{n) for real values of n in the interval [0,50] for the Ising model, that is the function n/2. 
Finally the solid red line gives the function n/(oo) for the sinh-Gordon model, that is a straight 
line passing through the origin which describes the asymptotic behaviour of the function n/(n) 
for n large. In the sinh-Gordon case, all functions have been computed for B = 0.5. (b) is a 
magnification of the lower left corner of the sinh-Gordon part of (a). 
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Figure 9: [Color online] The function n{f{0,n) 
various values of B. 



f{oo)) for integer n in the interval [2,50], for 



In principle we would only need to perform the integral in ()C.6p to obtain f{n) for any values 
of n. However this is highly non-trivial due to the complexity of the minimal form factor, as a 
function of n. We choose therefore to proceed in a different and more instructive way: we will 
instead find the natural analytic continuation of the function n/(0, n) numerically, as a large n 
expansion in powers of 1/n. 

First of all, it is an interesting exercise to try to determine analytically the precise slope 
of the line n/(oo) in the sinh-Gordon model. This is a relatively tedious but straightforward 
computation which we present in appendix[Dj The main result is the value of /(oo) as a function 
of B, which is given by 



/(oo) 



8192 



Bf{2 + Bf 



+ 



7r2 (4 
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in terms of the generalized hyper geometric functions defined in ()D.7|) and the variables 

5 5 



a 



3 
2 



5 



4' 



3 5 
2+4- 



(5.24) 



(5.25) 



As mentioned before, the data in Fig. [8] (b) also exhibit a clear deviation from the linear 
behaviour for small values of n, which is the region we are most interested in. Fig. [9] shows the 
values of the function n/(0,n) for integer n in the interval [2,50] with the linear part, that is 
n/(oo) subtracted, for various values of 5. By subtracting the linear part, the behaviour of the 
function for large n becomes more clear and appears to be dominated by a term proportional 
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to 1/re. This dependence can be made more precise by numerically fitting the various functions 
plotted above to a function of the generic form 



. ^ Oo ai 02 as 

Putin) = h^ + -^ + ^. 

n n' 



(5.26) 



The fact that only odd powers of n appear is a numerical observation. The table below contains 
the values of the constants ao, ai, 02 and as obtained numerically for various values of S, as well 
as the exact value of /(oo) from equation (j5.24p : 



B 


/(oo) 


ao 


ai 


02 


as 


0.02 


0.0952 


0.67(1) 


-0.17(2) 


-0.06(1) 


-0.02(3) 


0.1 


0.0972 


0.68(5) 


-0.20(5) 


-0.04(8) 


-0.02(3) 


0.2 


0.0994 


0.70(3) 


-0.24(2) 


-0.02(9) 


-0.02(7) 


0.3 


0.1013 


0.71(8) 


-0.27(7) 


-0.00(9) 


-0.03(1) 


0.4 


0.1030 


0.73(1) 


-0.30(8) 


0.01(0) 


-0.03(8) 


0.5 


0.1044 


0.74(2) 


-0.33(5) 


0.02(9) 


-0.04(4) 


0.6 


0.1055 


0.75(1) 


-0.35(8) 


0.04(5) 


-0.05(1) 


0.7 


0.1064 


0.75(8) 


-0.37(6) 


0.05(9) 


-0.05(6) 


0.8 


0.1070 


0.76(3) 


-0.38(9) 


0.06(9) 


-0.06(1) 


0.9 


0.1074 


0.76(6) 


-0.39(7) 


0.07(5) 


-0.06(4) 


1 


0.1075 


0.76(7) 


-0.39(9) 


0.07(7) 


-0.06(5) 



As an example, Fig. [THl shows the fit (15.260 for B = 0.5 together with the corresponding 
values of n{f{Q,n) — /(oo)) for integer values of n in the interval [2,50]. The fit (|5.26|) is in 
fact extremely good for all values of B. More precisely, we have checked that the values of 
n(/(0, n) — /(oo)) —p^t{n) for integer n are always smaller than 10~^ for all values of B included 
in the previous table. We are now in the position to compute /(I) numerically. Its value for a 
certain B corresponds simply to adding up all numbers in the corresponding row of the table 
above. The outcome of this computation is: 



B 


0.02 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


1 


/(I) 


0.5(1) 


0.5(1) 


0.5(0) 


0.5(0) 


0.5(0) 


0.5(0) 


0.4(9) 


0.4(9) 


0.4(9) 


0.4(9) 


0.4(9) 



As expected, all values obtained agree with the predicted value of 1/2 within our numerical 
precision. Hence, the entanglement entropy can be computed to 



SAirm) = log(em) + ?7sinh-Gordon - ^i^o(2rm) + O (e"^^'") . (5.27) 

We do not know how to evaluate the constant U yet, but the methods developed by 
S. Lukyanov [42j may be helpful. 



6 Generalization to theories with several particles and bound 
states 

So far we have been dealing with theories with a single particle spectrum and no bound states. 
In this section we wish to investigate how our results for the entropy can be generalized to 
situations in which the spectrum of the initial integrable QFT consists of £ > 1 particles and 
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Figure 10: [Color online] For B = 0.5, p{it{n) (|5.26p . and n(/(0, n) — f{oo)) at integer values of 
n, in the interval [2,50]. 

bound states are still absent. As before we will consider n copies of the theory and therefore it 
will be natural to label particles by two indices 

{a,i) with a = !,...,£ and i = l,...,n, (6-1) 

with the identification (a, i) = (a, i + n) for alH = 1, . . . , n and a = 1, . . . , i. Denoting by ^(a,i) 
some fundamental field of the theory related to particle {a,i), the exchange relations with the 
fields T and T can be written similarly as before: 

= 'r(x)^'(,,,+i)(y) x^>y\ (6.2) 

^'(„,,)(y)r(x) = T(x)^'(„,,)(y) < y\ (6.3) 

^'(„,,)(y)r(x) = t(x)^'(,,,_i)(y) x^>y\ (6.4) 

^'(„,,)(y)r(x) = t(x)^'(„,,)(y) x^<y\ (6.5) 

Denoting by Sa/3{d) with a,/3 = the two-particle S'-matrix of the original theory, the 

S'-matrix of the n-sheeted theory can be written as: 

S(a,i)if3,i){&) = Saf^{e) V a,P,i, (6.6) 
S{a,i)(/3,j)(.0) = 1 V a,/3,i,j with i j . (6.7) 

The form factors axioms (13. 5p for the operator T still hold, with /u the double index {a,i), 
jj, = {a,i + 1) and /2 = {a,i). If particles a,/3 in the original theory fuse to produce a bound 
state 7, then the 5-matrix has a pole on the imaginary line of the physical sheet, say iu^p with 
uj^p G (0,7r). Correspondingly, the form factors will possess extra poles with the requirements 

= r2X+7'^'""'""'(^i'---'^n-i) (6-8) 
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where the so-called three-point coupling is 

i^pY = -\ - Kp)So.m ■ (6.9) 

As usual, for finding solutions to these equations we must first construct minimal solutions 
of the two-particle form factor equations. They satisfy the equations 

-^min - ^min 1-^^' n)b(^oc,j){l3,k){'^) - ^min l^^^ " ^' (.O-^^J 

for all values of j, A;, a and (3. Prom the equations above it follows: 

^r|(M(Am)(^^„) = Flt'^^^^'^^'Ho^n) V (6.11) 

Flt''^^^'^\e,n) = F2^'^)(«'^)(27r(i-l)z-0,n) V a,/?,jV 1. (6.12) 

So, as before computing the minimal form factors of particles in the first sheet is sufficient to 
determine all minimal form factors of the theory. This minimal form factor must have no poles 
in the extended strip lm(0) G [0, 27rn] and satisfies a similar type of equations as for the case 
with ri = 1, 

^r|("'^)(/^'^)(^,n) = FT''^^-^'\-e,n)S^m = F^f ^^("'^^(-^ + 2.m,n). (6.13) 

From arguments completely analogous to those developed for the one particle case, provided 
that the S-matrix of the original theory admits an integral representation of the form 



Sap{0) = exp 



' dt , , . , fte\ 

I 75«M*)smh^-J 



(6.14) 



where QafiiG) is a function which depends of the theory under consideration with the property 
dafii't) = 9pa{t) (that is, parity invariance), the minimal form factor is given by 



(6.15) 



We are now in the position to obtain the full two-particle form factors by including appro- 
priate poles. In contrast to the single particle spectrum case not all two-particle form factors 

will have poles. The kinematic and bound-state residue equations ensure that the only singular- 
ities occur for form factors of the type F2^^'^'^^^"'''\6,n) for any values of i,j,a, or of the type 
j^-j for any value of a,P for which there is a non-zero three-point coupling T^^. If 
^aa = 0> then the pole structure of these form factors is analogous to the one found for the one 
particle case 

.ri(«,i)(.,i)... ^ {T)sm{--) F^ir''^^^^'\e,n) 

2nsinh(^)sinh(i^)i.2-^)(--i) (,,,„)■ 

In the presence of bound states, more factors need to be multiplied in order to account for the 
bound-state poles: 

pTliaMM^O^ = (^ + ^cosh(|)+Ccosh(^))F:r^(^-^)(g,n) 

(sinh {^) sinh sinh (%^) sinh (%t!) 

where A,B,C are constants that depend on a, /?, n. For all other form factors we have 

^r\{a,i){M^0^^pr\{a,imi)^^^^ ^^^^ r^;3 = 0. (6.18) 
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6.1 Computation of the entropy 



We can now proceed to the computation of the entropy along the hnes described in the previous 
section. The two-point function in the two-particle approximation is 

oo oo \ 
1 + ^^ / y d^l^^^2/a,/3(^12,n)e-^('"'^™^^^i+"^''™^^^2) (6.19) 
a,/3=l_oo -oo / 

where nia is the mass associated to particle type a. The new function fa,f3i(^,n) is given by 



r\{a,l)(p,i) 



i=l 



(6.20) 



,r|(a,l){/3,l) 



71-1 



{9,n: 



For general integrable QFTs it is also possible that the twist fields admit non-vanishing form 
factors both for odd and even particle numbers. If that would be the case there would be 
an extra contribution to the two-point function (|6.19p coming from the 1-particle form factor. 
However this would have no contribution to the entropy, since the 1-particle form factors have 
no singularities in 6 and therefore the 1-particle form factor contribution will exactly vanish at 
n = 1. 

For the computation of the entropy we have exploited the pole structure of the form factors. 
In fact, from our analysis it follows that only form factors containing poles will contribute to 
the final value of /a,/3(l), where /o,/3(n) is the natural analytic continuation of /a,/3(0,n) from 
n = 2, 3, . . . to n G [1, cxd). More precisely, our analysis showed that only when poles collide do 
we have a non-zero contribution to the analytically-continued sum above at n = 1 (recall that 
the first term vanishes at n = 1). Since < u^^ < vr, only kinematic poles can collide at n = 1, 
and this only occurs when ^ = 0. Selecting out those terms that contain kinematic poles, we 
obtain the constraint that /3 = a in (j6.19p for n ~ 1, which allows us to perform one of the 
integrals as before: 



(T(r)t(O)) « (T)2 I 1 + -^^ I d9fa,a{S,n)Ko{2rmcosli{e/2)) 



a=l 



(n ~ 1) (6.21) 



for fa,p{(^i'n) the natural analytic continuation of fa,i3{G,n). From here, the results of the 
previous section easily generalise to 



/a,a(l) 



d_ 

dn 



vr 



m 



n=l 



and 



SA{rm) 



1 ^ 

-\ log(em) + U--Y, Ko{2rm^) + O (e-^^-i) 



(6.22) 



(6.23) 



where nii is the smallest mass and we admitted the possibility of non-zero three-particle form 
factors. In particular, the leading exponential correction is given by the mass of the lightest 
particles, and the coefficient is proportional to the number of such lightest particles in the 
spectrum. 



27 



7 Conclusion 



In this paper we have developed a program which, taking full advantage of quantum integrable 
model techniques, allows for the computation of the entanglement entropy of a connected region 
of a quantum one-dimensional integrable system with respect to the remaining part of the 
system. The starting point of our computation is the well-known "replica-trick" . Our program 
relies in the realization that starting with a quantum integrable model and constructing a new 
model consisting of n non-interacting copies of the original theory, a new local quantum field 
theory is obtained which naturally possesses Z„ symmetry. Associated to this symmetry two 
twist fields T and T exist, whose two-point function is directly related to the entropy of the 
system. Through the "replica-trick" mentioned above, the entropy is the derivative with respect 
to n of the two-point function of the twist fields, evaluated at n = 1. 

Since the fields T and T are local fields of the n-copy theory, their two-point function can 
be computed by exploiting the form factor program for integrable models. More precisely, it 
can be expressed as a sum for different particle numbers over products of the form factors of the 
two fields involved. This gives a large-distance expansion: computing the two-point function in 
the two-particle approximation gives the leading behaviour at large distances. This expansion 
is in fact expected to converge rapidly, and this leading behaviour is often enough to describe 
the two-point function up to relatively small distances. 

The behaviour of the entropy as a function of the distance was already well-known for 
very short and very large distances. For r <C {m being the mass of the lightest particle) 
that behaviour is determined by the underlying CFT which describes the integrable QFT in 
the ultraviolet limit. Thus, the entropy can be computed explicitly by using CFT techniques 
|18t I19j . At large separations r ^ (in the infrared limit) the entropy is known to saturate 
to a constant value. The main result of our work has been to evaluate the first correction to the 
entropy in the infrared limit, providing therefore a description of the behaviour of the entropy 
in the intermediate region of values of rm. This correction is obtained from the two-particle 
contributions to the form factor expansion. In fact, the two-particle approximation provides 
both the saturation value of the entropy, coming from the disconnected part of the correlation 
function (the square of the vacuum expectation value of the twist fields), and the exact first 
corrections up to 0{e~^'^^), coming from the two-particle contributions. We also computed the 
exact value of the saturation in the Ising model and showed it to be in good agreement with 
previous numerical results [TTj . 

The most surprising result of our analysis has been to establish that the leading correction 
to the entropy at large rm is in fact a universal quantity, that is, it does not depend on the 
particular scattering matrix of the model we started with, but only on the spectrum of masses 
of the particles of the original theory. It is quite remarkable that the entropy should encode so 
explicitly crucial information about the theory both in its UV regime (the central charge) and 
in the IR regime (the number of light particles). 

We have deduced this result from general arguments and checked it explicitly for the Ising 
and sinh-Gordon models. The mathematical reason for this "universal behaviour" is clear, as 
the result is directly related to the presence of kinematic poles in the two-particle form factors. 
Only form factors having such poles do contribute to the final result for the entropy, and their 
individual contributions turn out to be theory-independent. The presence of bound state poles 
does not change our conclusions. However, we do not yet have a physical understanding of this 
result. 

It would be interesting to compute higher order corrections to the entropy coming from 
states with higher particle numbers. In addition, an interesting application of our main result 
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emerges, namely given a lattice model whose underlying QFT is unknown, a computation of 
the entropy could reveal the number of light particles of the theory by extracting the coefficient 
of the first exponential correction at large distances. Another interesting route is to compute 
the entanglement entropy of a disconnected region. This would involve higher-point correlation 
functions of the twist fields introduced here, and would be useful in verifying certain general 
properties of the entanglement entropy which follow from its interpretation as the number of 
"links" between the regions considered. 

An additional result of our work has been to develop the form factor program for branch- 
point twist fields. As a consequence of the particular exchange relations between these fields and 
the fundamental fields of the theory, the form factor consistency equations for branch-point twist 
fields are different from those associated to standard local fields. In particular, the crossing and 
kinematic residue equations are modified. Since here we have only been concerned with the two- 
particle form factors of branch-point twist fields, an interesting open problem remains, namely 
to find closed solutions to the form factor equations for arbitrary or at least higher particle 
numbers. From the form factor equations we expect these solutions to be given in terms of 
elementary symmetric polynomials of the variables e^'/". Also, it would be interesting to extend 
our analysis to non-diagonal theories, such as the sine-Gordon model. 

In the Ising model, it is possible, for even n, to have an explicit representation of the branch- 
point twist fields using their relation to U{1) twist fields in the free Dirac fermion model as shown 
in appendixlHJ This provides a way to explicitly evaluate Ising form factors, where many-particle 
form factors are obtained from two-particle form factors by the usual Wick's theorem. It also 
gives explicit representations of two-point functions in terms of the known Painleve transcendant 
representations of the two-point functions of U{1) twist fields. It would be interesting to study 
the Ising and Dirac models further from this viewpoint. 

Finally, the notion of twist field as introduced in this work is new in the context of integrable 
QFT and can still be generalized further. It would be very interesting, for instance, to study 
twist fields in integrable models which possess non-abelian internal symmetries (such as su{2)). 
Developing the form factor program for such fields is a possible future line of research. 

Acknowledgments: B.D. acknowledges support from an EPSRC post-doctoral fellowship, 
grant number GR/S91086/01. J.C. acknowledges partial support from EPSRC grant EP/D050952/1. 

A Form factors of branch-point twist fields from angular quan- 
tisation in the sinh- Gordon model 

It is possible to describe massive integrable quantum field theory using angular quantisation 
[64:\ I35j . Angular quantisation is obtained when the model is quantized on radial half-lines 
emanating from some fixed point in space. That is, the Hilbert space is some subspace of 
the space of field configurations on those radial half- lines, and time evolution is a rotation 
around that point. Let us consider the sinh-Gordon model, whose scattering matrix is given by 
(j3.25p . The angular quantisation of the model is obtained from free bosonic modes satisfying 
the canonical commutation relations 





with 



2 sinh ^ sinh 



n{2-B)u 
4 



(A.2) 



cosh ^ 
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along with the zero modes 

[P,q] = -i. (A.3) 

The hamiltonian of the system, generating rotation around the origin, is written in terms of the 
modes as 

K = / dv—-\.^K . (A.4) 

One can easily calculate 

[K, K] = -vK ■ (A.5) 

The angular Hilbert space, a Fock space over these modes, will be denoted by ^ = (BpJ^p- 
It is naturally decomposed into subspaces J^p with fixed eigenvalue p of the operator P (this 
determines the way the fundamental field behave when it approaches the origin 0, which is, in 
angular quantisation, the far negative infinity of space). 

The main advantage of this construction is that in some integrable models, an explicit 
embedding is known of the states in the usual quantisation scheme (with Hilbert space 7i) 
into the space of operators acting on J^. This embedding goes as follows [35]. One first identifies 
vacuum expectation values in TC with traces on J^, as a simple consequence of the change of 
quantisation scheme: 

TV (e-2^^6(A)) 
Tr (e" 

where on the left-hand sides A stands for an operator that is the representation on TC of products 
of local fields, and on the right-hand side l(A) stands for the representation on the angular- 
quantisation Hilbert space J- of the same fields. Then, there are operators Z{0) acting on 
such that products of them correspond to asymptotic states of TC: 

(vac|A|^i, ...,ek) = {{iiA)Ziei) ■ ■ ■ ZiOk))) . (A.7) 

The operators Z{0) are defined as follows. We have first 

^n^e)=:e~'^^'"'^"'"'"'-'^'"' : (A.8) 

for ?7 = lb, where the normal-ordering is with respect to the trace in (1A.6|) . and we define 

Z{e) = -iC [e'^A+(6' + iTT/2) - e'^^K'^e - i7r/2)l (A.9) 



(vac|A|vac) = ^ {{i{A))) (A.6) 



with 



C = — ■j^=^= exp 



sm^^ 



°° dt sinh ^ sinh 
t sinh t cosh | 



(A.IO) 



The representation l{A) on T of any field at the origin that is local with respect to the funda- 
mental field commutes with these operators; this constraint is expected to be sufficient to fix 
the set of all such fields in angular quantisation (which excludes twist fields - they are discussed 
below). 

For instance, form factors of exponential fields e*^"^ at the origin are obtained by choosing 
i{A) to be a projector on P = a, times the vacuum expectation value (e"*^). In particular, the 
identity operator is the projector on P = 0. Evaluating the traces is simple using 

{{KK')) = Y^^Ji^ + (A.ll) 



30 



obtained from cyclic properties of the trace. From the usual formula for evaluating averages of 
normal-ordered exponentials of free modes, 



» 



J dvdv'a{u)b{u'){{X^\,)) 



(A.12) 



we obtain 



exp 



^i+r/if JA^^ (^2 + ry2f ])) = 

dt sinh ^ sinh / Oi- 

-'^ViV2 / T r— 771 cosht 1 + z 

Jo t smhrcosh^ \ n 



2(^1 -'72) 



(A.13) 



Two-particle form factors of exponential fields are then given by 



(vac|e"'^|ei,e2> = (e'"^)^' 



'?l,2=± 



'i + ^i?)A^^ 



ITT 



^2+r/2y ))> (A.14) 



and the constant C ensures that the value of the kinematic residue is correct. 

Angular quantisation is also useful for studying twist fields. Since the "space" (or equal-time 
slices) of angular quantisation is just the half line, a twist field associated to a symmetry action 
cr is just the operator for the symmetry action itself on the angular Hilbert space: 



i{%{0)) = {%)a^ 



(A.15) 



and this also projects on the subspace with P = 0. Note that the symmetry action applied 
on a field, [ajr,i[0{x ^ 0))], inside the trace on the angular Hilbert space, does produce the 
symmetry transformation of the field 0{x): since there is a commutator, there are two branches 
inserted at slightly different angular times, and the two branches cancel each other except around 
the point x. Also, these two branches can be deformed into one branch from x = —00 to a; = 00, 
so this is really the same object as the operator for the symmetry action on TC. When there is 
no commutator, that is when just cjjf is inserted inside the trace, then this does not produce 
symmetry transformation, rather it corresponds to the insertion of a local field, the twist field. 

We now consider the n-copy sinh-Gordon model: the new angular quantisation Hilbert space 
is just (B^=iJ'^"'\ where all J^^"''^ are isomorphic to J-', and the asymptotic state operators are 
just Zj{6), for j = 1, 2, . . . , n. They are built as before out of bosonic modes Xj^u, j = 1, 2, . . . , n 
that commute for different values of j. The local twist field T = %j associated to the symmetry 
c • J ^ J + 1 mod n is then just given by (jA.lSh . 

Form factors of T are then 



(vac|T|^i 



with 



/Mi,...,Atfe 



{T)c'^{{zM...z^(ek))). 



(A.16) 



(A.17) 



We now take the operators Z^{9) to be composed of exponential operators normal-ordered with 
respect to this new trace, in order for the computation of the trace of their products to go 
as before. Such a change of normal-ordering, for exponential of free modes, just changes the 
normalisation, and the operators A^(^) all get the same normalisation change. The constant C 
has been introduced in order to account for this, and can be fixed by the requirement that the 
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correct value of the kinematic residue is obtained from the two-particle form factors evaluated 
below (it could also be determined directly by a Bogoliubov transformation). 

Using the fact that the symmetry acts like aj^Xj^^a^^ = Aj+i^i, (with j + n = j) and cyclic 
properties of the trace, we can evaluate the following averages: 

-2{j-l)7VI^ f( \ 

{{X,,uXiy)). = '^('^ + (j = 1, 2, . . . , n) (A.18) 

Note that with j = 1, the average is the one obtained by changing the angle around the origin 
to 27rn. This is naturally expected: the branch points represented by these branch-point twist 
fields are just negative-curvature conical singularities with an angle of 27rn. For j > 1, the 
formula says that in evaluating form factors with two particles belonging to different copies, 
one can replace Xj^u by e~'^^^~^^'^'^ Xi^^- From formula (lA.Sp . it is clear that this is equivalent 
to changing the particle I, of type ji, to type 1, and shifting its rapidity 9i to 9i + 27ri(j; — 1), 
as long as the resulting integral representation for the form factor is convergent after this shift. 
This is so if the shift keeps the rapidity in the strip lm(0i —62) G [0, 27m], so that this shift is 
allowed if it is the first particle, I = 1 with rapidity 61, that is changed. Hence, we find 

Fp\ei - 62) = Fp'\9i -02 + 2Tri{j - 1)) (A.19) 

which is in agreement with p. lip (this equation holds for the whole form factor, not just the 
minimal form factor). 

Hence it is sufficient to evaluate form factors involving particles on the same copy. The 
evaluation of traces involved in calculating these form factors goes as above, and we find 



/•°° dt sinh ^ sinh ^2-^ ( q^-Q^ \ 

-2??i% / -T r-— ri — cosht n + i -{r]^ - ri^) 

Jq t smh nt cosh ^ \ ir Z 



(A.20) 



exp 
which gives 

(vac|T|^i, ^2)1,1 = {T){CCf ^-^^''^''^^^((A?^(^i + ^if)A?(^2 + r?2f))). 

»7l,2=± 

= i^)iCCf (-2F:S.\0) + I + I ] (A.21) 

Till 

(with 9 = 01- 02) where F^^^\0) is given by (13.280 . This can be seen to reproduce ()3.18p if we 
choose 

-iC./sin2^ 

iCC? = , . . , (A.22) 

2 (cos 2^ + sm 2^ — 1 j 



using the properties 
1 



, , . = /(2v^^n - 9)f'X{0) (A.23) 



with 

^ sin^(l+..)sin^(2+u;) ^ ' ^ 

and 

w = n H . (A.25) 

TT 



32 



B Vacuum expectation values of branch-point twist fields in the 
Ising model 

In this appendix, we use the relation between branch-point twist fields in the n-copy Ising 
model and the n independent U{1) twist fields in the n-copy free massive Dirac theory, in order 
to deduce vacuum expectation values of the former from the known formulas for those of the 
latter. Note that this derivation is very similar to that employed in |21j . 

The symmetry a : i i + 1 mod n associated to the branch-point twist field in the n- 
copy Ising model can certainly be diagonalised on the particle eigenstates. It would be natural 
to associate the resulting diagonal elements with elements of a U{1) symmetry group, and to 
interpret the branch-point twist field as a product of the corresponding n independent ^7(1) 
twist fields. This, of course, is unnatural in the n-copy Ising model, since there is no U{1) 
symmetry in the individual copies. The trick is to further double the n-copy Ising model in 
order to make it into an n-copy Dirac theory. Denoting the fundamental real Majorana fermion 
fields by ipa,j^ '^bj-i ^a,ji ^b,j for j = 1,2,... ,n and fundamental Dirac spinor fermion field 

= ( ^.^'^ 1 , we have the relations 



Each copy of the Dirac fermion has a U (1) symmetry: it is the symmetry under rotation between 
copies a and h of the Ising model, which occurs because the Ising model is quadratic. In the n- 
copy Dirac theory there is also an extra symmetry under SU (n) transformations of the multiplet 
*i \ 

• • • I , again because the theory is quadratic. Diagonalising the branch-point twist field Toirac 

^-n / 

in the Dirac theory, which is done with a SU (n) transformation as is seen below, it is possible to 
interpret the new basis as n new independent Dirac fermions. It is then possible that the branch- 
point twist field can be simply written as a product of n independent U{\) twist fields acting 
on these independent fermions. This is useful for the Ising model, because by uniqueness of the 
branch-point twist field (characterized by its basic branch-point property and its dimension), 
we have (recall that the central charge of the Dirac theory is 1) 

Toirac =%(^Tb (B.2) 

where 7^ and 7^ are the branch-point twist fields in the copies a and b of the n-copy Ising model 
respectively. 

In order to go into the details, we need to be more precise about the SU{n) transformation 
that diagonalises Toirac- It would be natural to expect that the Fourier transform "^^'^^ = 
^ 'Y^=i e'^^^^'^^j is the appropriate transformation, but there is an important subtlety: the 
various copies of the Dirac fermions commute amongst each other, whereas to form new fermions 
with linear combinations of them, they need to anti- commute. This is clear from our choice of 
scattering matrix Sjk{0) = 1 amongst different copies j ^ k. Of course, from the viewpoint 
of the Hilbert space this is only a choice of basis, and it is possible to define a new basis with 
scattering matrix —1 amongst different copies. There are many ways of doing that, all with the 
same result. We will choose the following (here with a two-particle example for simplicity): 

\o /) \ac _ / I^'l^2)jij2 Jl < h /T^ o\ 
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Written using annihilation and creation operators for the basis with upper index ac, fermions of 
different copies anti-commute (and this is the only change). In the new basis, the symmetry a 
acts as follows: 

^ac ^ f ^i+i j = 1, . . . , n - 1 , . 

^ \ -^f j = n . ^ ^ 

The extra minus sign is understood as follows: in a non-zero correlation function, there are an 
even number of fermion fields. Let us order them in increasing copy label. There will be, say, 
k fields with copy label n at the right. Applying the symmetry makes them into copy 1, and 
they are the only fields that break the order of increasing copy label. Bringing them back to the 
beginning gives a minus sign if they are of odd number. This is cancelled by the extra minus 
sign above. 

Writing a in matrix form, X]fc=i ^jk'^k = '^^j; we find that the eigenvalues A, Yl^=i ^jk^k = 
Xvj, satisfy 

A" = -1 ^ A = e"^ , p odd . (B.5) 
The explicit SU{n) transformation is 

n 

q,(p) — ^ ^g-ii-pi/n _ g-i7rp{l+i/")^ _ (B.6) 
i=i 

Hence, we may expect that, for n even, 

77. Tl_ 

rD.ac=no£^^^noits^^ (B.7) 

q=l q=l 

where Oa^ is the U{1) twist field acting non-trivially only on the fermion fields 'J'^^^, associated 
with the U{1) element e^'^*", for a G [0,1]. The dimension of these twist fields is a^, and the 
choice of the odd values of p above forming the product is dictated by the requirement of having 
the lowest total dimension. The dimension of the product of fields on the right hand side of 
(|B77l) is 

n 

which agrees with the CFT prediction (|2.8p for central charge c = 1. By uniqueness of the 
branch-point twist field, the factorisation above must be the correct one. We checked in the case 
n = 2 that our form factors agree with formula ()B.7p along with the known form factors form 
C/(l)-twist fields [43]. For n odd, similar calculations show that there are many lowest-dimension 
factorized operators, and that their dimension does not agree with the CFT prediction (it is 
higher). Hence, we expect that in this case, no such factorisation exists. This is not a problem, 
as we are interested in the analytic expression for the vacuum expectation value, which can be 
obtained from even values of re. 

We can now read off the expectation values of Oa from 



V /m\"^ f °° dt f sinh^ at n n.\ /rre\"^ 1 ,^ , 

<''"> = (y) T(iw?r-°= ) = (y) f^-') 

which gives (15.140 (here G{z) is Barnes' G-function). 
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C The function f{e, n) 

We now derive the full analytic continuation f{6, n) in an integral representation. For simplicity, 

7~l 11 

we will assume that F2 {0,n) is at = and that it vanishes exponentially as \6\ 00, as 
is observed in the Ising and the sinh-Gordon models. Let us write 

n-l 

f{e,n) = Y,sie,j) (c.i) 

"Till ( 7"l 11 \ * 

with s(6l,j) = F2 ' (-61 + 27rij,n) F2 ' {-e-2mj,n). Consider now the following contour 



integral: 

dz 

7rcot(Trz)s(9,z) (C.2) 

/c ^vrz 

where C is the closed rectangular contour with vertices (n — iL, n + iL, iL, —iL). Then, s{6, z) 
decays exponentially as Imz itcxD, and the contributions from the horizontal segments vanish 
as L — > cxD. Thanks to the quasi-periodicity of the integrand, s{9,z + n) = S{9 — 2Triz)S{9 + 
2TTiz)s{9, z), the contributions from the vertical pieces amount to 

dz 

/ {S{9-2TTiz)S{9 + 2mz)-l)TTCot{T^z)s{z) . (C.3) 

J-ioo 27r? 



The residues from the simple poles a.t z = j , j = 1, . . . ,n — 1 sum to f{9,n). There are also 

and z = n-^±^ 

2 2m 2 2m 



simple poles of s{z) at z = i ± 2^ and z = n — ^±^1- the kinematic singularities. They 



evaluate to (for real 9) 

tanh(^Q lm(^Fp^\-29 + i7r,n)-F[^^\-2e + 2mn-i7r,n)^ . (C.4) 

At 6 = this gives —1/2 using the kinematic residue equation, as it should. For any non-zero 
6, it vanishes as n ^ 1, like (n — 1)^. As an analytic function of 6, there are simple poles at 
6 = ibz7r(n — 1), which are those that give the main contribution as n — > 1 for 9 near to 0. These 
poles are 

icot(^) icot{^) 
2{9 + i7r{n - I)) ' ~ 2(6* - i7r(n - 1)) 

which indeed gives the behavior (|5.5p . with /(I) = 1/2. 
The full analytic continuation can now be written: 

fi9,n) = tanh(^Q Im (^Fp^^{-29 + iTT,n) - F^^^ {-29 + 2Trin - iTT,n)^ (C.5) 

{S{e-f3)-S{e + P))Fp^\-9 + P,n) (^Fp^\e + P,n)y d(3 . 

In particular, at 9 = 0, this specializes to 

/(n) = \-^ r Im(5(-/3))coth \FP'\l3,n)\^dp . (C.6) 
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D Computation of /(oo) for the sinh-Gordon model 

In this appendix we present a detailed derivation of the equation of the hne that provides the 
large n behaviour of the function /(O, n) (and by analytic continuation that of f{n)) in the 
sinh-Gordon model. The computation can be carried out as follows: the limits of the various 
factors entering the two-particle form factor at ^ = 2'Ki[j — 1) are given by 



lim Fjj^ {2m{j -l),n) 

n~*oo 



sm 

lim 

2n sinh i '-^^^^ sinh ( '-^^^^ 



2n 



[J 



^' r(j-f)r(i 
i r(i-|)r(i + |) 
2r(i-f)r(i + f)' 



1 + 1) 



7r(2j-l)(2j-3)^ 



1 , 

2 4 y 



for j <C n, 



for j <C n, 



(D.l) 



where for the first two functions, the limit can be easily evaluated by employing the integral 
representation ()3.28p of the minimal form factor. Putting all factors together we obtain 



lim Fp'^^{2TTi{j - 1), 

I— ►oo 

4(j - 1) r 



nj 



|)r(j-i + f)r(i-|)r(i + |) 



vr(2,-i)(2,-3) rO--|)r(j-i + f)r(i-|)r(i + |) 



(D.2) 



for j <C n. In order to compute /(O, oo) we still need to perform the sum over j of the function 
above. Since the sum goes up to j = n and in (|D.2p we have assumed that j <C n, it is 
not guaranteed that (|D.2p will provide a good approximation for the whole range of values 
j. It turns out that (jD.2p is the right function to consider for two reasons: first, computing 

(2'7rz(j — 1), n) numerically for various values of j and n one quickly realizes that it is only 
non-zero for values of j around 2 or around n; second, due to the periodicity of the form factor 
at ^ = 0, it turns out that the value of the function (ID.2P does not change when j ^ n — j + 2. 
Therefore, if we call the function ()D.2p k{j), we can replace Yl^=2 ^0) ~^ ^ Sj=2 ^(•?) n. — > oo. 
We then obtain 



/(0,oo) 



4 ; 



r 1 + 



r 1 



4)T 



MB), 



with 



A{B) 



oo 



(J - 1) 



l)2(2j-3)2 



r(j-i-f)r(,-i + f) 



2 4 ; 



Employing the definition of the Pochhammer symbol 



(a) 



In 



T(a + n) , , , , , 
= V/ N = a « + 1 « + 2 . . . a + n - 1 , 
Pa) 



(D.3) 



(D.4) 



(D.5) 



we can rewrite the sum above as 



AiB) 



f)r(i + f 



B-] 



r(2-|)r(i + |) 



i=2 



(2j-l)2(2i-3)2 



4^i-2 



1 + 



(D.6) 
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and carry it out analytically. The result is given in terms of generalized hypergeometric functions, 
which are defined as 



We find 



F 



oi, a2, . . . , Up] z 
61, 62, •••,&(? 



(Qi)fc(a2)fc • ■ ■ (Qp)fc z'' 
^0 {bi)k{b2)k...{bg)k k\ 



(D.7) 



with 



E 

i=2 



(i - 1)^ 



(2j-l)2(2j-3)2 



(2 4 )j-2 (2 + 4 /j-2. 



9 



+ 



75 



a/3 



1225 



a/3 



2 / 



|,|,2,a + l,a + l,/3 + l,/3 + l;l 



2' 2' ' 



(8 + ^)^(^-10)^ 
(10 + 5)2(5-12)2 



1,«; + 1,(7 + l,cj + 1 
7-^6 



2 , 2 ) 1) CK) CK) /5) /5) 1 
2 ' 2 ' 



5 5 
2' 
9 

2' 2' ' 



2, 2,3,a + 2,a + 2,/3 + 2,/3 + 2;l 
^ ^ K + 2,Av + 2,cj + 2,o- + 2 



a 



2 4' ^ 
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